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MODELING  THE  BIDIRECTIONAL  REFLECTANCE 
AND  SPECTRAL  ALBEDO  OF  A  CONIFER  FOREST* 


Crystal  Barker  Schaaf  Alan  H.  Strahler 

Geophysics  Directorate,  Phillips  Lab,  Center  for  Remote  Sensing,  Boston  University, 

29  Randolph  Rd,  Hanscom  AFB,  MA  01731,  USA  725  Commonwealth  Ave,  Boston,  MA  02215,  USA 

ABSTRACT 

Analysis  of  aircraft  imagery  and  ground  measurements  acquired  for  a  spruce  forest  stand 
in  Howland,  Maine  as  part  of  the  1990  Forest  Ecosystem  Dynamics  Multisensor  Aircraft 
Campaign  (FEDMAQ,  show  that  the  Li-Strahlcr  geometric-optical  forest  canopy  reflectance 
model  can  estimate  spectral  bidirectional  reflectance  factors  and  spectral  albedo  with  some 
success.  In  this  model,  the  vegetated  surface  is  treated  as  an  assemblage  of  partially  il¬ 
luminated  tree  crowns  of  ellipsoidal  shape,  and  through  geometric-optics  and  Boolean  set 
theory,  the  proportion  of  sunlit  and  and  shadowed  canopy  and  background  is  modeled  as  a 
function  of  view  angle.  The  model  is  driven  by  ground  measurements  of  spectral  reflectance 
and  tree  crown  shape,  sire  and  spacing.  Atmospherically  corrected  multiangular  radiance 
measurements  from  the  Advanced  Solid  State  Array  Spectroradiometer  (ASAS)  fit  the  shape 
of  the  modeled  reflectance  function  quite  well.  Integration  of  *he  reflectance  function  yields 
a  spectral  albedo,  which,  when  extended  to  the  full  solar  spectrun,  agrees  with  pyranometer 
measurements  obtained  at  the  site. 

1.0  INTRODUCTION 

The  anisotropy  exhibited  by  forest  canopies  is  determined  by  the  canopy  structure,  shape  and  density, 
shadowing  patterns,  multiple  scattering  interactions  within  the  canopy,  leaf  angle  distributions  and  the 
specularity  of  the  foliage.  Although  dense,  uniform  canopies  have  primarily  been  modeled  with  homo¬ 
geneous  plane-parallel  or  three-dimensional  radiative  transfer  models  (Veihoef,  1984;  Gerstl  and  Simmer, 
1986;  Kimes  etal.,  1987;  Kimes,  1984;  1991;  Nilson  and  Kuusk,  1989;  Myneni  etal.,  1989;  Verstraete  et 
al.,  1990),  simpler  geometric-optical  models  have  shown  some  success  in  modeling  non-uniform  and  par¬ 
ticularly  sparse  forests  (Li  and  Strahler,  1985;  1986;  1992a;  b;  Strahler  and  Jupp,  1990a;  b;  Abuelgasim 
and  Strahler,  1993).  Because  sparse  canopies,  or  forests  with  significant  height  variations  produce  com¬ 
plex  shadowing  patterns,  the  interplay  between  the  sunlit  and  shadowed  components  of  the  canopy  is  the 
primary  determiner  of  the  radiation  reflected. 

The  Li-Strahler  geometric-optical  model  exploits  these  shadowing  effects  by  modeling  a  scene  or  pix¬ 
el  as  an  assemblage  of  ellipsoidal  tree  crowns.  Geometric-optics  and  Boolean  set  theory  are  used  to 
determine  the  areal  proportions  of  shadowed  and  sunlit  canopy  and  shadowed  and  sunlit  background  as¬ 
sociated  with  a  solar  zenith  angle.  Independent  characteristic  signatures  for  each  component  are  weighted 
by  these  areal  proportions  and  used  to  determine  the  spectral  directional  reflectance  factor  of  the  canopy 
at  any  viewing  angle.  Although  these  component  signatures  do  not  include  an  explicit  treatment  of  diffuse 
irradiance,  canopy  multiple  scattering  or  leaf  specularity  and  are  applied  uniformly  to  the  areal  propor¬ 
tions,  they  are  reasonable  approximations  for  clear  atmospheres  and  low  to  medium  solar  zenith  angles. 
The  geometric-optical  model  does  incorporate  the  effects  of  mutual  shadowing  or  the  obscuring  of  tree 
crowns  by  one  another  and  takes  into  account  the  impact  bright  unshadowed  tree  tips  have  on  the  canopy 
reflectance  at  very  high  view  angles.  By  performing  a  hemispherical  integration,  the  model  also  provides 
an  instantaneous  hemispherical  reflectance  (or  spectral  surface  albedo). 

♦Presented  at  the  25th  International  Symposium,  Remote  Sensing  and  Global  Environmental  Change, 
Graz,  Austria,  4-8  April  1993. 
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The  availability  of  directional  radiances  from  the  Advanced  Solid  State  Array  Spectroradiometer 
(ASAS)  obtained  during  the  1990  Forest  Ecosystem  Dynamics  Multisensor  Aircraft  Campaign  (FED- 
MAC)  in  Howland,  Maine  makes  a  validation  of  the  model  possible.  Atmospherically-corrected  ASAS 
reflectances  are  typically  obtained  in  at  least  the  principal  plane  and  cross-principal  plane  of  the  sun. 
They  can  be  compared  (both  in  shape  and  magnitude)  to  modeled  bidirectional  reflectance  factors  (BRFs) 
when  the  model  has  been  carefully  calibrated  with  accurate  tree  shape  measurements  and  component  sig¬ 
natures  measured  in  the  field  at  the  time  of  the  ASAS  overpass.  Although  sky  radiance,  multiple  scatter¬ 
ing  and  specular  effects  are  not  explicitly  incorporated  in  the  Li-Strahler  model,  they  are  implicit  in  the 
field-measured  component  signatures. 

In  addition  to  validating  the  model  BRFs,  the  spectral  hemispherical  reflectances  computed  by  the 
geometric-optical  model  can  be  combined  (Brest  and  Goward,  1987;  Brest,1987)  to  produce  a  surface  al¬ 
bedo  that  can  be  validated  with  pyranometcr  measurements  obtained  on  the  site  at  the  overpass  times. 

2.0  PROCEDURE 


2.1  ASAS  IMAGES 

At  9:12  EDST,  8  Sept  1990,  the  ASAS  instrument  was  flown  on  a  NASA  C-130  aircraft  at  4603  me¬ 
ters  altitude  in  the  principal  plane  of  the  sun  over  the  FEDMAC  spruce  site  .  This  site  is  in  the  vicinity  of 
an  instrumented  meteorological  tower  maintained  by  the  University  of  Maine  at  Orono  (45°  21.21  ’N,  68° 
44.49’ W).  ASAS  images  are  obtained  at  seven  look  angles  from  +45  °  to  -45°.  The  red  band  15  (band 
center  644.6  nm)  and  near-infrared  band  24  (band  center  773.5  nm)  were  selected  for  this  study.  A  fairly 
uniform  region  of  vegetation  was  selected  around  the  tower  and  the  mean  brightness  for  this  area  was  ob¬ 
tained  from  each  image.  This  brightness  was  then  transformed  to  radiance  with  the  NASA-provided  ra¬ 
diometric  resolution  factors.  Finally  the  values  were  atmospherically  corrected  using  the  Liang  and 
Strahler  (1992)  model  as  parameterized  by  aerosol  optical  depths  measured  on  the  site  at  the  time  of  the 
overpass.  The  cross-principal  plane  images  were  obtained  somewhat  later  at  9:40  EDST  and  subjected  to 
the  same  processing. 

2.2  FOREST  CHARACTERISTICS 

In  1989  investigators  from  Boston  University  surveyed  the  spruce  stand  near  the  tower,  and  used  vari¬ 
able  radius  plot  sampling  (Dilworth,  1977)  to  obtain  tree  height,  diamctcr-at-breast  height,  crown  width, 
and  height-to-crown  distance.  These  data  were  supplemented  with  tree  heights  collected  by  University  of 
Virginia  researchers.  Based  on  the  sample  of  40  trees,  the  parameters  used  as  input  to  the  geometric- 
optical  model  were  basal-  area-weighted  means  of  hcight-to-ccntcr-of-crown  (9.9  ±  2.4m),  crown  radius 
(2.2m),  crown  vertical  radius  (3.6m),  and  density  (1161  trces/hectarc). 


In  1990,  at  the  time  of  the  ASAS  overpass,  radiances  were  collected  using  a  handheld  radiometer 
(Spcctron  Engineering  SE590)  for  use  as  tire  sunlit  and  shadowed  component  signatures  required  by  the 
model.  Bands  97  (643.7  nm)  and  141  (774.1  nm)  were  used  to  match  the  band  centers  of  the  ASAS  chan¬ 
nels.  Sunlit  and  shadowed  crown  measurements  were  taken  within  half  an  hour  of  the  overpass  time. 
Spcctralon  panel  measurements  were  also  taken  so  that  reflectance  factors  could  be  computed.  The  sha¬ 
dowed  and  sunlit  background  signatures  were  composites  of  radiances  from  litter,  fems,  shrubs,  moss  and 
grasses.  Unfortunately,  the  background  signatures  were  measured  about  2  hours  after  the  overpass.  How¬ 
ever,  given  the  fairly  large  solar  zenith  angle  (58.79  °)  and  the  high  density  and  crown  coverage  (84%), 
the  crown  components  would  be  the  primary  contributors  to  the  scene  reflectance  rather  than  the  back¬ 
ground  components.  So,  the  later  background  signatures  were  used  (with  no  attempt  to  correct  for  solar 
zenith  angle).  A  three-component  version  of  the  model  was  used,  so  that  the  shadowed  component  is  the 
average  of  the  background  and  canopy  signatures.  The  final  component  signatures  used  in  the  model  for 
the  red  spectral  region  were  a  sunlit  canopy  of  0.0349,  a  sunlit  background  of  0.0503,  and  shadows  of 
0.0017.  For  the  NIR  spectral  region  values  used  were  a  sunlit  canopy  of  0.3969,  a  sunlit  background  of 
0.3248,  and  shadows  of  0.0375. 
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2.3  ALBEDOS 


Upward  and  downward  pointing  Eppley  PSP  pyranometers  were  positioned  at  a  height  of  25.8  meters 
on  a  boom  extending  2.4m  west  of  the  meteorological  tower  at  the  spruce  site.  These  data  can  be  radoed 
to  produce  instantaneous  lull-spectrum  surface  albedos  (.3-3pm).  During  the  times  of  the  ASAS  over¬ 
passes,  these  surface  albedos  range  from  0.1 17  to  0.1 1 1. 

3.0  RESULTS 

3.1  MODEL  VERSUS  ASAS  REFLECTANCES 

Using  the  forest  characteristics,  the  component  signatures,  and  the  appropriate  solar  zenith  angle,  the 
geometric-optical  model  was  used  to  simulate  the  spruce  site  at  the  time  of  the  ASAS  overpass.  The  red 
and  NIR  bidirectional  reflectance  distribution  functions  (BRDFs)  are  displayed  in  Figure  1.  The  three  di¬ 
mensional  BRDFs  are  displayed  in  a  rectangular  coordinate  system  where  each  view  angle  in  the  hemi¬ 
sphere  is  taken  as  a  pair  of  polar  coordinates  and  transformed  onto  the  x-y  plane  as  a  vector  of  unit  length. 
The  corresponding  reflectances  are  then  plotted  along  the  z  axis.  The  hotspot  occurs  where  the  viewing 
angle  and  solar  zenith  angle  coincide  thus  concealing  shadows.  The  bowl  occurs  when  the  view  angle 
moves  opposite  the  solar  zenith  angle  and  more  and  more  of  the  shadows  are  revealed.  The  width  of  the 
hotspot  is  governed  by  the  elliptical  shape  of  the  crowns.  The  shape  of  the  bowl  is  governed  by  the  densi¬ 
ty  of  the  forest  and  the  amount  of  mutual  shadowing  that  is  taking  place.  By  extracting  results  along  the 
principal  plane,  the  modeled  near-infrared  BRFs  can  be  compared  to  the  ASAS  BRFs  (Figure  2).  The 
model  predicts  brighter  values  in  the  forward  scattering  directions  than  were  detected  by  the  ASAS  im¬ 
agery.  However,  the  general  shape  and  magnitude  were  quite  good.  The  modeling  was  repeated  with  the 
same  component  signatures  but  with  the  solar  zenith  angle  that  corresponded  to  the  cross-principal  plane 
flight  (54.54°).  In  the  cross-principal  plane  direction  (Figure  3),  the  model  overestimates  the  impact  of 
the  mutual  shadowing  by  creating  more  of  a  valley  shape  than  the  shallow  bowl  shape  revealed  by  the 
ASAS  imagery.  However,  the  magnitude  is  quite  reasonable. 

3.2  MODEL  ALBEDO  VERSUS  MEASURED  ALBEDO 

An  integration  over  the  viewing  hemispheres  of  Figure  1  result  in  a  spectral  hemispherical  reflectance 
or  surface  albedo.  Brest  and  Goward  (1987)  suggest  that  the  red  spectral  model  albedo  is  representative 
of  a  0.526  proportion,  that  the  near-infrared  model  albedo  is  representative  of  a  0.362  proportion  and  the 
mid-infrared  is  representative  of  a  0.1 12  proportion  of  solar  radiation  incident  at  the  surface.  Therefore, 
red,  near-infrared,  and  mid-infrared  spectral  albedos  can  be  weighted  by  these  proportions  and  added  to 
estimate  the  broadband  value.  Since  Uiere  were  no  mid-infrared  radiometer  measurements  made  to  use  as 
component  signatures,  no  mid-infrared  model  runs  could  be  made.  However,  Brest  and  Goward  (1987) 
use  0.5  of  the  near-infrared  value  as  an  appropriate  approximation  of  the  mid-infrared  value.  Using  these 
assumptions,  full  spectrum  model  surface  albedos  can  be  computed  and  compared  with  the  pyranometer 
measurements  taken  at  the  site.  The  geometric-optical  model  produced  values  of  0.1025  and  0.1102  at 
the  times  of  the  ASAS  overpasses,  while  the  pyranometer  measured  values  during  that  time  range  from 
0.1 17  to  0.1 1 1  (Figure  4).  These  model  values  were  quite  close  in  magnitude  to  the  measured  values. 

4.0  CONCLUSIONS 


The  Li-Strahler  geometric -optical  model  was  validated  with  data  collected  from  a  spruce  forest  in 
Maine  during  the  1990  FEDMAC  data  collection  effort.  Although  modeled  near-infrared  BRFs  agreed 
quite  well  with  ASAS  reflectances  along  the  principal  plane  in  both  shape  and  magnitude,  a  comparison 
of  the  cross-principal  plane  values  reveal  that  the  model  tends  to  overestimate  the  reflectance  at  mid-to- 
high  viewing  angles  by  slightly  overemphasizing  the  impact  of  mutual  shadowing.  Model  spectral  al¬ 
bedos  were  combined  to  produce  full  spectrum  albedos  that  could  be  compared  with  pyranometer  meas¬ 
urements  taken  at  the  site  at  the  same  time  as  the  ASAS  overpasses.  Given  the  simplicity  of  the  model 


and  the  difficulties  in  obtaining  realistic  and  timely  component  signatures,  the  close  agreement  found 
between  the  measurements  and  model  results  is  quite  encouraging. 

Such  success  demonstrates  the  adaptability  of  the  geometric-optical  model  but  also  highlights  its 
dependence  on  accurate  and  timely  sunlit  and  shadowed  component  signatures.  Work  is  underway  to 
produce  a  hybrid  model  that  treats  canopy  scattering  as  a  function  of  path  length  and  gap  probability  (Li 
and  Strahler,  1988;  Strahler  and  Jupp,  1990)  and  combines  the  scattering  with  the  geometric-optical  ef¬ 
fects  of  the  trees  within  the  canopy.  Such  efforts  will  result  in  a  model  that  obviates  component  signa¬ 
tures  and  does  not  suffer  from  a  sensitivity  to  site  and  time  specific  spectral  measurements. 
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Figure  1.  Red  and  NIR  geometric -optical  model  bidirectional  reflectance  distribution  functions  at 
the  time  of  the  ASAS  principal  plane  overpass  (SZN  58.79  °,  GMT  13:12) 
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Figure  4.  Measured  (-)  and  Modeled  (♦)  Surface  Albedos  at  ASAS  Overpass  Times. 
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